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We derive a factorization theorem for the differential distributions of electroweak 
gauge bosons in Drell-Yan processes, valid at low transverse momentum, using the 
Soft-Collinear Effective Theory (SCET). We present the next-to-leading logarithmic 
(NLL) transverse momentum distribution for the Z-boson and find good agreement 
with Tevatron data collected by the CDF and DO collaborations. We also give predic- 
tions for the Higgs boson differential distributions at NLL based on a factorization 
theorem derived in earlier work. We derive formulae for all quantities needed to 
study low transverse momentum production of color-neutral particles within the ef- 

[- — , fective theory. This effective field theory approach is free of Landau poles and can 

be formulated entirely in momentum space. Consequently, our results are free of the 

i— i Landau-pole prescriptions necessary in the standard approach. 
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I. INTRODUCTION 

The study of low transverse momentum (jpx) production of electroweak gauge and Higgs 
bosons plays an important role in many studies within and beyond the Standard Model (SM). 
The measurement of the PF-boson mass at hadron colliders requires an understanding of W 
production at low pt pp. The search for a Higgs boson decaying into two W bosons requires 
a jet veto to remove the large ti background [2]. This cut effectively restricts the Higgs boson 
to the low-px region. In the kinematic limit of low transverse momentum, large logarithms 
of the form ln(M/pr), where M denotes the mass of J*,W, Z, h, spoil the perturbative 
expansion based on the strong coupling constant as- The logarithms must be resummed to 
all orders to obtain an accurate prediction. This resummation has been extensively studied 
in the literature [5HTS] . The standard approach utilizes a Fourier transform from momentum 
space to impact-parameter space to decouple emissions of multiple gluons while maintaining 



momentum conservation [3]. This introduces a Landau pole arising from evaluating the 
strong coupling a s (l/b\) for large impact parameters, b± —¥ oo. The Landau pole must be 
dealt with for any transverse momentum, even for px 3> ^qcd- The resummed exponent in 
6j_-space also does not vanish when px becomes large, which potentially introduces numerical 
instabilities in the matching to fixed-order QCD at high transverse momenta. An alternative 
approach which addresses these limitations is worth pursuing. 

In a previous paper we performed an analysis of transverse-momentum resummation 
for the example of Higgs-boson [19] production using the Soft-Collinear Effective Theory 
(SCET) [20H22] . We derived a factorization theorem for the differential distributions of 
the Higgs boson. The resummation of large logarithms in this approach is performed using 
renormalization-group (RG) evolution in the effective theory. We introduced several new 
functions which encode the emission of soft and collinear gluons within the effective theory, 
which we labeled impact-parameter beam functions (iBFs) and inverse soft functions (iSFs). 
Our SCET formulation has several advantages over the standard approach. The matching 
from full QCD onto SCET PT at the hard scale M followed by running to lower scales is done 
in momentum-space, not in impact-parameter space as in the standard method. This has 
the effect of stopping the RG evolution at the scale jjlx ~ Pt, the natural scale of SCET PT 
describing the dynamics of soft and collinear px emissions. The factorization theorem can 
be written entirely in momentum space which is equivalent to the statement that the in- 
tegral over b± can be performed analytically without running into the Landau pole. As 
a result, for perturbative values of px, the transverse momentum distribution is predicted 
entirely in terms of perturbative functions and the standard initial-state PDFs. Only for 
non-perturbative values of Pt, one obtains a new non-perturbative function in SCET PT which 
is field theoretically well defined and has a computable anomalous dimension. This is in con- 
trast to the standard approach where a prescription is introduced to avoid the Landau pole 
even for perturbative values of px- Since an analytic integration over impact parameters is 
possible, we avoid instabilities that occur in the standard approach in the matching of the 
resummed exponent to the fixed order result needed in the region of high px- 

Our goal in this manuscript is to extend our previous work in several ways. We formulate 
the factorization theorem for low px production in SCET PT to account for electroweak gauge 
boson production in addition to Higgs production. We present one-loop expressions for all 
iBFs and iSFs needed to study resummation to the next-to-leading logarithmic (NLL) order 
of accuracy. The two-loop results for the iBFs and the iSF are required for a complete 
resummation to next-to-next-to-leading logarithmic (NNLL) accuracy. We compare the 
SCET formulation to the standard Collins-Soper-Sterman (CSS) [7J approach to transverse- 
momentum resummation. We present numerical results for the Higgs and Z-boson px dis- 
tributions at NLL. We compare the effective theory predictions for Z production with the 
experimental measurements at the Tevatron (231 123] , and find very good agreement with the 
data. This demonstrates that our SCET-based approach provides a promising alternate way 



of studying transverse momentum distributions at hadron colliders. We outline future direc- 
tions for the study of transverse momentum resummation within the effective field theory 
framework. 

Our paper is organized as follows. In Section [II] we pedagogically review the factorization 
theorem derived in our previous work [TH] and present its extension to electroweak gauge 
boson production. The details of this extension are presented in appendix |Aj All analytic 
results for the matching coefficients, iBFs and iSFs required for phenomenology to NLL and 



partial NNLL accuracy are presented in Section III The structure of the RG running in 



the effective theory, which resums large logarithms of the form In (M/pr), is discussed in 



Section IV Simple analytic expressions for the resummed cross sections valid through NLL 



are shown in Section [Vl We discuss the relationship between the various quantities appearing 



in the SCET approach with those appearing in the CSS formulation in section VI, and show 
the consistency of the methods through NLL. We discuss what further work must be done 
to establish the relationship to higher orders. Numerical results for Higgs production and Z 



boson production are shown in Section [VII and the agreement with the Tevatron data for 



Z production is demonstrated. Finally, we conclude in Section VIII 



II. REVIEW OF THE FACTORIZATION THEOREM 

We begin by summarizing the content and derivation of our previously- studied factoriza- 
tion theorem [TH], and present its extension to the case of electroweak gauge boson produc- 
tion. The details of this extension are presented in appendix |XJ The derivation and result 
of our factorization analysis are shown schematically below: 

d 2 a 



| (match QCD to SCET PT ) 

PS\C®(O SC et)\ 2 

\, (SCET soft-collinear decoupling) 

~ H ®B n ®B n ®S 

\, (zero-bin and soft subtraction equivalence) 

~ H ^Bn^Bn^S' 1 

\, (match SCET PT to SCETa qcd ) 

~ H ® [In ® Xn <8> S' 1 ] ®fi ® fj. 

v v ' 

s 

In the first stage of the analysis, full QCD is matched onto an effective field theory 
which contains fields with the following momentum scalings: 

p n ~ M(t] 2 ,1,t]), p n ~ M(l,r] 2 ,r]), p s ~ M(r),r],ri), 7] ~ -^, 



corresponding to the n-collinear, n-collinear, and soft modes respectively. M denotes 
the mass of either an electroweak gauge boson or the Higgs. We consider only the 
leading operators in the rj expansion, which encode the most singular emissions coming 
from soft and collinear particles. The relevant operators for Higgs and electroweak 
gauge boson production are respectively 

0(ui,u} 2 ) = g^hT {Tr S , «(0.B£ L ) wl S£Sft(0.B£ I ) wa S& }, 
OfKa*) = &W)n,„ 2 T[S n r« i S} l ](W^ i ) n , u;i , (2) 

where the T^ 1 denotes the Dirac structure, the index K runs over the vector and 
axial-vector Dirac structures, and the index i runs over the quark flavors. The B^ ± 
and B£ ± fields denote collinear-gluon field strengths [25] dressed with collinear Wilson 
lines in the n and n directions respectively, £j denotes a collinear quark field of flavor 
i, and W denotes a collinear Wilson line. The ui^ are label momenta that give the 
large light-cone components of the collinear fields. 

• Using the soft-collinear decoupling property of SCETpr, the matrix element in 
SCETpx is decoupled into the n and n collinear iBFs, B n and B n respectively, and a 
soft function as seen in the fifth line of Eq. (JT1). 

• The B n and B n iBFs are defined with a zero-bin subtraction |26j. Using the equivalence 
between zero-bin and soft subtractions, explicitly shown at one-loop in our previous 
work [TJJ] and studied elsewhere in the literature [2"TH2"5] . we can write B n B n S — )■ 
BnBnS^ 1 where the iBFs B Ujfi are defined without the soft zero-bin subtraction and 
S- 1 is the iSF. 

• In the last step, for p? 3> Aqcd-, the iBFs are matched onto standard PDFs via the 
schematic matching equation 

B n ,n = -Ln,n ® J ■ (<Jj 

The matching coefficients X n ^ are grouped together with the iSF to form a Transverse 
Momentum Function (TMF) Q as shown in the last line of Eq. (fTj). The TMF encodes 
the physics of the soft and collinear emissions from initial-state partons in the effective 
theory. For p T ~ Aq C d, the OPE in Aq C d/pt breaks down so that the perturbative 
matching in Eq.pl) is no longer valid. However, in this case one can view Eq. (pi) simply 
as an equation that defines new non-perturbative functions X nfi . Thus, for pt ~ A-qcd, 
Q appears as a new non-pertubative TMF with a well-defined field-theoretic definition 
and computable anomalous dimension and can be modeled and extracted from data. 

The above derivation relies on the cancellation of Glauber mode contributions, as in the 
standard approach [7], to the final observable which measures the pt of the final-state color- 
neutral particle. An explicit demonstration of this cancellation using an effective field theory 
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language remains to be shown. For Higgs production, our previous work [19] showed that 
the differential distribution in transverse momentum and rapidity of the Higgs is given by 

d 2 Oh n 2 f 1 dx\ f 1 dx 2 f 1 dx\ f 1 dx 2 

dp 2 T dY = 4(iV c 2 - 1)2Q2 J Q ~^J ^ J xi lc[ ] X2 ^ 

x H h (xi, x 2 , n Q \ ht)Qh(x u x f 1 ,x 2 , x' 2 ,p T , Y, n T )fi/p{x'n /J>T)fj/p(x' 2 , fir), (4) 

where the we have introduced Q% , theTMF for Higgs boson production, which is given by 

//■ r a2i. 

dt$ / dt n / — — J (b ± pr) g^gpu 

X Zn;gi(zj> t n > &1> /J>t) Z™gj(—T> ^> b±, f*r) 

x 1 x 2 

n r\ —Y / 2 2 n 

-—,x 2 Q-e Jv T + m 2 h n- 

x 2 Q v xiQ 



x S gg 1 (x 1 Q-e Y Jp' T + m 2 h --^ } x 2 Q-e Y \/v T + m\ - -^, b±,fPr). 



(5) 



Detailed definitions of the various objects appearing in the above equation are given in 
sections |IHC and HID The integrations over impact parameter bj_ can be explicitly per- 



formed to produce an expression written completely in momentum space [19] . Similarly, the 
integrations over the residual light-cone momentum components t^,t^ can be performed. 
For clarity, the explicit expression that results after performing the indicated integrations is 
shown in Section [V] 

In this paper, we also give a detailed derivation of electroweak gauge boson transverse 
momentum and rapidity distributions in appendix |A} The final result is given by 

d 2 a 7T 2 f 1 , f 1 , f 1 dx' f 1 dx' 

" dx\ I dx 2 



dp T dY N c J J J X1 x 1 J X2 x 2 

x H q z {x x x 2 Q 2 , fi Q ; /Mr) G qrs {x x , x 2 , x'^x'^pr, Y, /tfr)/ r (x'i, /j, T )f s (x 2 , /z T ), 

(6) 

where the TMF function Q qrs 

g qrs (x 1 ,x 2 ,x' 1 ,x' 2 ,p T ,Y,fi T ) = / -pr^ J o[b±Pr] dt^dt^X n . qr (^,t^,b ± ,fi T )X n . gs (^-,t^,b ± ,fi T ) 



x S qq \ Xl Q - e Y sj V 2 T + M 2 - ^,x 2 Q - e~ Y \Jp 2 T + M 2 - ^,&±,/*r) 



and detailed definitions of the various objects can be found in appendix [X} One can straight- 
forwardly generalize these results to distributions that are differential in the final state lep- 
tons. 

We briefly address the appropriate choices for the scales \iq and fix that appear in these 
equations, as this issue arises later when detailed numerical results are presented. These 



(7) 



scales should be chosen to minimize logarithms that appear in the expressions for the hard 
function H and the TMF. A detailed study of this issue was given in our previous work [T9] ; 
we summarize the results of this analysis here. The scale \xq which appears in the hard 
function should be chosen as \xq ~ X\X 2 Q ~ Mz- The scale \xt appearing in the TMF 
should be chosen to be p? ~ px- We note that the TMF is independent of the quantities £+, 
t^ and b± which are integrated over. The TMF is only sensitive to the transverse momentum 
and rapidity constraint on the final state Higgs or electroweak gauge boson making px ~ Pt 
the natural scale choice. 

III. FIXED-ORDER RESULTS 

In the following sections we derive fixed-order results for the various components of the 
factorization theorems in Eqs. Q and (J6J). We provide mostly outlines of the necessary 
perturbative calculations, as details of the derivations were already presented in our previous 
work [Tj5] for the case of Higgs production. 

A. Matching Coefficients 

We begin with the matching coefficients that arise when matching the full QCD current 
onto operators in SCET PT . The vector-boson current in full QCD is 

■fc;/. = E( J & + J &)' (8) 

i 

where Jj. 4 and J^ 1 are vector and axial- vector currents for the i-th quark flavor 

Jz& = 9vQil^Qi, J£l = #a&7m75&, (9) 

and g l v and g\ are the vector and axial vector couplings of the i-th quark to the vector boson 
being studied. We have explicitly used the Z current in these equations, but the results hold 
identically for 7* and W bosons as well. These currents are matched onto effective operators 
in SCET PT as 

4i = j dun j du 2 C K > ji (u u w 2 , riO* 1 ^, w 2 , /i), (10) 

where the index K takes on the values K = {V, A} and the indices i, j run over the quark fla- 
vors 1 . C K '^ % is the Wilson coefficient of this matching. The effective operators O^ 1 (l0i, uj-i, p) 
are given by 

Of ( Wl ,w 2 ,/i)(x) = {iw)n^T\s\£fs n \(yv^ n ^ (11) 



1 There is also a pure gluon SCET PT operator that should appear on the RHS of Eq.(10). However, the 



contribution of this operator vanishes [3D] for Drell-Yan processes and has thus been left out. 



with the Dirac structures T^ written as 



^vi 



9vX 



i-i- ■ 






(12) 



S nt n denote soft Wilson lines 

r-0 



S n (x) = Pexp 



n ■ A s (x + sn) 



bn{X) 



P exp 



W 



n ■ AAx + sn) 



(13) 



Note that in the Dirac structures T^ q in Eq.(12) only the perpendicular components of 
the index /i contribute due to the equations of motion of the collinear fields in SCET PT : 
^£n — 0, ijl^n = 0. At tree level, the Wilson coefficient arising from matching the current 
onto the SCET PT operator is just unity: 

C F;ji(0) (w b w 2 ,/i) = C*t i W(u u u 2 ,n) = 6 ij . 
The result of the one-loop matching is given by 



(14) 



a s C F 



ln 2 



P 



-UJ1U2 — It 



3 In 



F 



-bJ\bJ2 — It 



TX 

~6~ 



(15) 



The expression for this Wilson coefficient up to order a 2 , can be found in Refs. 
The hard Wilson coefficient H%(ui, o>2, H>Q\ ^t) = Hz(u>iU>2, /J>Q] A*t) that appears in the 
factorization theorem in Eq.d6| is given by Eqs.(A20), (A21), and (A27). We refer the 
reader to [19] for a similar analysis for the hard coefficient H h (xiX 2 Q 2 , (jlq; Ht) that appears 
in the factorization theorem of Eq.Q for Higgs production. 



B. Quark iBF 



In this section we give the results for the one-loop calculations of the quark iBFs that 
appear in the factorization theorem for electroweak gauge boson production. The special 
case of the quark iBFs with the transverse coordinate b± = was computed at one loop in 
[M] . Our computation closely follows the techniques we established in our previous work [TJJ] . 
We compute the n-collinear and the n-collinear quark iBFs by inserting a complete set of 
states \X n ) and \X n ) respectively as 

' I r -"■.'■" > lnAf....\\ JIT. In )\ \.,)( X J _ 



B*(x,t,b±,[i) 



2xn ■ pi 
1 



<!l> e^ b ~ ^(pi\^W n {b-,b x )\X n )(X n \h(n ■ V - x n ■ Pl )W^ nq (0)\ Pl ) , 



2x n ■ p2 J 47r 



Air 



sCn 



- 6 2Q: 



b+ 






spin 



rfi 



(W^ q )n(y)\X n )(X n \5(-u - Pfi)(Z q W)*(0)\ \p 2 ), 

(16) 



e 



•oooqoooo' 



i 

FIG. 1: Example diagrams for the real emission of a single parton contributing to the next-to- 
leading order quark iBF. The purple cross denotes the collinear Wilson lines associated with the 
Xn held. We note that the momentum p\ is incoming on the left-hand side of the cut and outgoing 
on the right. 



and then computing the product of matrix elements. The tree level and and virtual cor- 
rections are obtained by a perturbative calculation with the choice of the vaccum state 
\X n ) — |0). The virtual corrections are all scaleless and vanish in pure dimensional regular- 
ization. The real emission contributions correspond to choosing final states \X n ) to contain 
one or more partons. Example diagrams for the real emission of a single parton are shown 
in Fig. [T] The matching of the iBFs onto the PDFs is given by 

- 1 dz 



B q n (x,t,b±,[i) 



ln;qqi(-,t,b ± ,{l) f q '(z, /l) + X n - qg (-, t, 6± , fi) fg(z, fi) 



(17) 



with an analogous equation for the n-collinear iBF. The PDFs are defined in the usual 
fashion as 

f q (z^) = \ J2 (Pi\L q W n (0)h(n-V-zn- Pl )W^ m {0)\Pi), 



fg(z,V) 



initial pols. 
zn ■ Pi 



d-2 



E (Pi\[Bi a Ab^b ± )6(V-zn. Pl )Bi ±a (0)}\ Pl ). (18) 



initial pols. 



In the rest of this section we give results for the n-collinear iBF. Analogous results hold for 
the n-collinear iBF. At tree level the n-collinear iBF is given by 



B~f\x,t,b±,fi) = S(t)S(l-x). 



(19) 



At the next order, the emission of a single parton into the final state from the iBF has two 
contributions: 



BfU(x,t,b±,») = K qR{1 \x,t,b ± ,ri + B« 9R V(x,t,b±,»), 



(20) 



\qqR{l) 



iqgR(X) 



where Bn and B n correspond to the first and second diagrams in Fig. [l| respectively. 
The contribution of a single gluon emission by an initial state quark to the n-collinear quark 
iBF is given by the first diagram of Fig. [TJ All other diagrams do not contribute if a physical 
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polarization sum is used for the final-state gluon. At the level of the integrand, the result 
for this diagram is 



B?? R U(x,t,b ± ,n) 



9 2 C F 



d d k 



t 



x 



2xQ J (2n) d 

k^,n v + K v n^ 



-9(k Q )5(k 2 )5[Q(l -x)- k-]5[— - k + ] 



Ak-b i 



xQ 



fiV 



n ■ k 



Tr 



jfi 



"(*-*)£(*- W 



(k- Pl y 

(21) 



where Cp = 4/3 and we have introduced Q = n ■ p\. This can be computed in pure 
dimensional regularization to give 



BfM(x,t,b±,l*) 



a s C F e^ / b\t(l-x) 

ro-^i i — e; — 



2tt r(l-e) 



Ax 



1 r/i 2 i [ --' 

7 



'1 - e)x € (l - x) 1 - 6 + 



2.v 



l+e 



'l-x) l+e 



(22) 



We now expand this expression in e and use the matching condition of Eq. (17). In dimen- 



sional regularization, the only contribution to the PDF is at leading-order, as higher-order 
corrections are scaleless. Since f(x) = 5(1 — x), we have 



Xn m (x, t, b±,/j,)= B%*(x, t, b±, //) 



- finite part in dim. reg. 



(23) 



We have also performed a calculation where collinear divergences are regulated by introduc- 
ing an off-shellness for the initial quark, and have obtained identical results for the matching 
coefficients. We find the following expression for the matching coefficient: 



2&(:M,&.L,/i) = ^r{m 



2n 



+ S(t) 



5(1 — x) hire + (1 — x) 

6 1 — x 



ln(l — x) 



/;/ 



/' 



1 — x 
2 ^ 1 + x 2 



'l + x 2 ) + 



f' z 



ln(t//i 2 



5(1 -x) 



t J + (l — x) + 



FA1 



b\t(l - x) 
Ax 



(24) 



Another contribution to the quark iBF comes from setting the final state \X n ) to a single 
quark, as shown in the rightmost diagram of Fig. [TJ This contribution matches to the gluon 
PDF, and generates the qg — > Vg partonic channel that contributes to electroweak gauge 
boson production. The integrand level expression for this diagram is given by 



B q n 9RW (x,t,b ± ,(,) 



d d k 



2 Qx j w-t^X* 2 )^ - *1W - ») - *-]^^r^ 



~9fiu + 



Q 



Tr 



'h-f<jh.-h)\v*-tirt 



(25) 
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In dimensional regularization this becomes 



B^\x,t,b ± ^) 



a, e 



£1E 



47iT(l-e) 

1 r/i 2 i 1+e 



Ml-e: 



Wt(l-x) 



Ax 



l' z 



„2+e/ 



2x' + \\ - x)- e - 2x l+e 1 - x)~ e + (1 - e)x £ (l - x) 



(26) 



Expanding this result in the limit that e — > and keeping only the finite remainder, we 
derive the matching coefficient 



T%\,.1.h A _,n = g^(t)[{(l-x) 2 + x 2 } In 1-^ + 1 



+ 



IL 

/jH t J 



{(1 - xf + x 2 } FJl 



b\t(l-x) 
Ax 



(27) 



Gluon iBF 



We now consider the calculation of the gluon iBF. Although parts of this computation 
were already presented in Ref. [19], we include them here for notational consistency and 
completeness. Note that in the following we use slightly different definitions for the iBFs 
compared to that in [19J. The n-collinear gluon iBF is given by 

db- 



B% afi (x,t,h ± ,n) 



4tt 



_g 2 xQ 



J2 J2^\[9 B t±0(b-M\X n ) 

initial pols. X n 



X 



(X n \6('P-x 1 n-p 1 )gBi ±Q (0)}\p 1 }, 



(28) 



where on the right-hand side we use the Fourier transform conjugate variable t/(xQ) com- 
pared to t/Q used in [19]. The contributing diagrams are shown in Fig. [21 In dimensional 
regularization and with a physical polarization sum used for the gluons, these are the only 
contributions. The matching equation for the gluon iBF onto the PDFs is given by 

B g nvx0 (x,t,b±,t*)= f ~ {^£M±,rif g (z)+^&t,b±,l*)f 9 {z)}. (29) 



The tree-level expression for the iBF is given by 

B$W{x,t,bi.,l*) = 9 2 gfS(t)8(l-x). 



(30) 



The virtual corrections to the iBF are scaleless and vanish in pure dimensional regularization. 
At the next order beyond tree-level, the real emission of a single parton from the iBF has 
two types of contributions: 

BS;\x,t,b ± ^) = B^\x,t,b ± ^) + B^ 1 \x,t,b ± ^), 



(31) 
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FIG. 2: Example diagrams contributing to the next-to-leading order gluon iBF. The purple cross 
denotes the collinear Wilson lines associated with the B n ± field. We note that the momentum p\ 
is incoming on the left-hand side of the cut and outgoing on the right. 



where -B 99 a i and B^ 9 a a correspond to the emission of a real gluon and a real quark in the 
final state respectively. We begin by considering the contribution arising from single-gluon 
emission into the final state, which corresponds to the leftmost diagram in Fig. [2] The result 
can be expanded in terms of two form factors, 



K- a p(x, t, bx,fi) = Ff'^x, t, b ± ,n)g'f + Jf e (x, t, b ±} //) 
where T\p. are given by 

gg;e „ 2 a sC A ^ 1 V 



a/3 



9? + (d ~ 2) 



b«b\ 



H 



9 



7v T(l-e)/j 2 [ t 
1-x 1 



l+e 



x 1+£ (l-x)- £ Fi 1-e; 



bjtjl - x) 

Ax 



x < 1 — x 



x- 



1 — X 



~C99'> e 
• r 1 



9 



2 a s C A e^ E 
Ati r(3-e) 



t 



1 — x 



,2-6 



X' 



o^i (3 -e; 



b 2 ± t(l - x) 

Ax 



(32) 



The quantity J-" 2 is finite, and we can immediately set e to zero. T\ must be expanded in 
distributions. We do so and drop the pole terms, as explained in the section on the quark 
iBF. For the gluon iBF we utilize the matching equation 






wiih't'bj-'M&v)' 



(33) 



which leads to the result 



X nL(x, t, b±,/j)= ^(x, t, b±, n)gf + J% 9 (x, t, b±, (j) 



a/3 



9T + (d-2 



Iff If, 



(34) 
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with 



Tl 9 = 9 2 Ca- ~S(t) 5(1 -x)+ 5(1 - x) 

71 VI 



+ 5(t)x 
1 



1 — X 
I - x 4 — — ] ln(l — x) 
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-Ix 



(35) 



There is another contribution to the gluon iBF coming from radiating a quark into the 
final state, which comes from the rightmost diagram of Fig. [2j This matches onto the quark 
PDF, and generates the matching coefficient X n . gq . As with B^ 9 a/3 , this can be expanded in 
two form factors, leading to the generic form 

K aP {x,tM,n) = J c ! q (x 1 t,b±,M + J c i q (xAb±^)[g^ + (d-2) b ^]. (36) 

o 

After performing the integrations and expanding in e as in the previous contributions, and 

rt 

bU(l-x) 



utilizing the matching in Eq. (33), we find the following expressions for the two form factors: 

;i-a;) 2 + i 
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m 
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2 oi s Cf 



2tt [h 



l r^ 
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)*i(i; 



Ax 
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'1 — x) 2 + 1, 1 — x 1 — X 
'- In 2 



x 



X 



X 



TOT/ + U \ 2 a sC F (l-x) 2 b 2 ± 



2n 



Ax 2 



oFi (3; 



b\t(l - x) 
Ax 



(37) 



D. Inverse Soft Functions 



We now discuss the computation of the inverse soft functions that appear in the fac- 
torization theorem for both Higgs and electroweak gauge boson production. They serve to 
remove an overcounting of soft emissions that occur when the iBFs are inserted into the 
factorization theorem, and are needed to correctly match the fixed-order cross section. The 
relevant iSF for gauge boson production is 

db + db- hQlh+e ^ b - s _, {h+x ^^ 



S qq (wi,u 2 ,b±,fi) 



16tt- 



(38) 



where S qq (b,fi) is a vacuum matrix element of soft Wilson lines given by 
S„(b,n) = ^Tt(0\f[S} l S n ](z)\X M )(X.\T[£liS n ](Q)\0). 



(39) 
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We have inserted a complete set of final states \X S ). The relevant formulae for the iSF in 
Higgs production are 

db + db- 

'"99 



Sj{(o 1 ,u 2 ,b±,ij) 



16tt- 



3 » + Oi/2 e *-Oa/25-i(6+ > 6- > 6_ L ) ) 



j-max j— jmax j.+ 

0—1(0)/^ 2^_ l _n l Ji_ 

qq l x 2 Q ' Xl Q " 



S gg (b,fi) = J2(°\ f [ Tr [SnT D StS n T c Slj (b)\ \X S )(X S \T [Tr {S n T c S^T S± j (0)J |0). 

(40) 
The tree level result for the quark iSF is given by 

&±,a0 = iv c xix 2 g 2 (y(c ,a, -*n)*(*r*-C). 

(41) 



where we have used the same arguments in the iSF that appear in the factorization theorem 
so that 

t% ax = x x x 2 Q 2 - x 2 (M 2 - u), t™ x = x x x 2 Q 2 - Xx(M 2 - t). (42) 

u and t are the hadronic Mandelstam invariants that appear in the factorization theorem 
derived in appendix [X| The result for the gluon iSF at tree level can be obtained via a simple 
scaling of the quark result: 



c-i(o) = N c~ 1 c-i(o) 

99 4JV C 11 ■ 



(43) 



As with the iBFs, the contributions from the higher-order virtual corrections with \X S ) — 
|0) are scaleless, and vanish in dimensional regularization. The only contributions come 
from real-emission diagrams. Examples are shown in Fig. [3] for both vector boson and Higgs 
production. For gauge boson production, the contribution of real gluon emission to the iSF 
takes the integrand level form 

d d k 6(A; 2 )0(£; o )e*--^ 



j-max j— j-max j-+ 

C-lR(l)/hi v n_ l _ri L 
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X 2 Q 
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xiQ 
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Evaluating this in dimensional regularization and dropping the pole terms yields 

^ 2 1 



S -1R(1) 
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FIG. 3: Example diagrams contributing to the next-to-leading order iSFs for both electroweak 
gauge boson production (left) and Higgs production (right). The lines at each vertex schematically 
denote the soft Wilson lines associated appearing in the definition of the iSFs. 

We have used Q 2 = X1X2Q 2 to simplify the notation. The result for the iSF in Higgs 
production at this order can be obtained via the scaling 

C-LR(l) = ( N C - l )°A c-ifl(i) (A6) 

99 4N C C F qq ' K ' 



IV. RUNNING 

In the section we briefly summarize the structure of the running of the various objects that 
appear in the factorization theorems for the transverse momentum and rapidity distributions. 
For a more detailed discussion of the running structure see Ref. [TH] where the case of Higgs 
production was studied in great detail. The overall structure of factorization is similar for 
both Drell-Yan and Higgs production and is schematically characterized by a hard function 
H, a transverse momentum function (TMF) Q, and the standard initial state PDFs. The 
TMF Q is a convolution over two iBFs and an iSF. It is only the specific forms of the hard 
and transverse momentum functions and the the type of parton PDFs that dominate that 
differ in Drell-Yan and the Higgs production processes. In this section we summarize the 
running structure for the case of the Z-boson distribution function. 

The evolution equations for H^(Q 2 ,fi) are diagonal in flavor (for a more detailed discus- 
sion see [201 EU ES] ) , and one can write 

^H% = 1 H% HI (47) 

where the anomalous dimension 7^ has the form 

lHl = n[a s ]hx% + ^ H [< (48) 

The first term proportional to In % is known as the cusp anomalous dimension and the 
second term is the non-cusp piece. The quantities r#[a s ] and 7#[a s ] have perturbative 
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expansions in a s of the form 



T q H [a s 



a s r Hq 
47T 



r h 4- 

1 + 



47T 



1 2 



,ffg 



r^ + 



7hM 



¥ T0 + 



47T 



1 2 



7i +••• 



(49) 



The matching from QCD onto SCET PT is the same as for the study of threshold resummation 
in SCET, and the results for the anomalous dimensions can be obtained from previous 
studies [ST|[5B]. For resummation to NLL accuracy, the following coefficients of the expansion 
are needed: 



-,-ffq 



ii 



Hq 



8Cf, 


Hq 

7o = 


-V2C F , 


8C f 


Y67 vr 2N 
A 9 3, 


) c* - t* 



(50) 



The solution to the evolution of the hard function between the scales \iq ~ Q and \x? ~ Pt 
and takes the form 



H q z (Q 2 , u.q] u. T ) = U H i (Q 2 , u,q, u, T )H g z (Q 2 , u. Q ), 



(51) 



where U H i (Q 2 ,[iq, /j,t) denotes the evolution factor and sums the logarithms of Q 2 /p T - 

The TMF function Q qrs that appears in the factorization theorem as seen in Eq.(p| is 
evaluated at the scale \xt ~ Pt- The PDFs are evolved via standard DGLAP equations up 
to the scale \xt ~ Pt, summing up the remaining logarithms of Aqcd/pt- For the running 
between fig ~ Q ~ M and /xy ~ Pt, one can also consider running of the iBFs and the iSF 
individually above the \xt ~ Pt scale. It was shown in [19] that the combined convolution 
running of the iBFs and the iSF cancels the running of the hard function as required by the 
scale invariance of the cross-section. The running of the iBF can be written as 



(JL—B*(x,t,b ± ,(i) = dt'^{t-t',fi)B q {x,t',b ± ,fi), 



(52) 



where at one loop, the anomalous dimension is given by 



7f (1) (^) 



a s (n)C F 
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f' Z 



t J 



+ \m 



which is the same as what was found for the iBF with b 
dimension of the iSF S^ 1 is determined by the equation 

1 
Q~ 2 



(53) 
in 1 341. The anomalous 
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(54) 



and at one loop is given by 
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(55) 
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V. ANALYTIC EXPRESSIONS FOR RESUMMED CROSS SECTIONS 



Using the fixed-order expressions and renormalization-group evolution of the hard func- 
tion derived in the previous sections, we can derive the explicit expressions for the differential 
cross sections which resum large logarithms for low transverse momenta. We begin by con- 
sidering the case of Higgs boson production. The relevant factorization formula is shown in 
Eq. Q. The running of the hard function Hh was described in Section IV For the trans- 



verse momentum function defined in Eq. ([5]), we must plug in the perturbative expansions 
for the relevant iBFs and iSF. These were derived in Section IIII1 We note that if tree-level 
expressions are used for both beam functions and the soft function, then the phase-space 
constraints force px = 0. Therefore, the NLO expressions for one of these functions must be 
used. 

We derive here the 0(a s ) expressions for the differential cross sections, which correspond 
to the leading-order result for the px spectrum. Upon plugging in the expressions for the 
iBFs and the iSF, the integrals over b±, £+, t^, x±, and x 2 can be performed. The differential 
cross section can be written in the general form 

2 ij 



dp T dY 



\ dx ± dx 2 f i /p(x 1 )f j /p(x 2 ) Q - ' , (56) 

. . J dudt 

h3 



where the subscripts h, Z refer to the Higgs and Z boson production respectively. We have 
introduced the usual partonic Mandelstam invariants Q, t, and u. There are three relevant 
partonic channels that contribute to Higgs production: gg — > gh, qg — ¥ qh, and qq — > gh. 
The partonic differential cross section for the gg initial state can be written as follows: 



d 2 a{ 9 7T 2 


(a s y\ Q 2 


dudt 192t 


> 2 V TX ) | tU 


tu 


u u 2 
. Q 2 Q 4 . 
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x 5(Q 2 


+ t + u-m 2 hi 
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Q 2 

U h (Q 2 ,VQ,HT) + jr 
tu 
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t t 

1 + ^ + ^ 

Q 2 Q 4 



U h {m h -u,fj J Q,^ T ) 



TT , 2 x , 2t 2 + 2m 2 + 3tu + 6Q 2 m 2 h 
U h {m h - t, n Q , Ht) + -7- 



(57) 



The four terms appearing in this result have a clear origin. The first term arises from <S g }, 
while the second and third come from the NLO expressions for the iBFs. The last term 
results from matching the expressions to the fixed-order QCD result. In the limit that the 
p T of the Higgs becomes large, and therefore that ut approaches uq, the evolution factors 
Uh — > 1, and the regular fixed-order QCD result is obtained. Of the remaining channels, the 
qg initial state has only a contribution from the NLO iBF, while the qq has no contribution 
from soft or collinear emissions at this order and comes entirely from matching to the fixed- 
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order QCD result. The explicit expressions are 
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(58) 



For vector boson production, we for simplicity explicitly write the result only for on-shell 
Z-boson production with the leptonic phase space integrated over. Results for W and 7* 
production can be obtained through simple modifications of this formula, as can the results 
when the leptons are treated differentially. In this case, two partonic channels contribute, qq 



and qg. The result for the cross section can be written exactly as for the Higgs in Eq. (56) 
The partonic channels can be written as follows: 
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+ 2 
Q 2 Q 2 \ 



6(Q 2 + t + u-M 2 z ), 



(59) 



with the vector and axial couplings as defined in Eq. (pL The remaining partonic channel 
erf? can be obtained from <r| 9 by the interchange of the u and t variables. At the order to 
which we are working, the explicit form of the evolution factor U H i (£, /xq, /iy) appearing in 



Eq. (59) is given by: 



U h « z {£,Hq,Ht) = |exp{25(/i Q ,/i T ) -a r (/iQ,/i T )ln(-^//iQ) -a 7 (/iQ,ii T )} 
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ri 
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Cq 
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8/3o 



lnr, 



(60) 



where T ^ = r o ^/2,7 l? = 7 9 /2, and r = ocs{^t) / o/.s{^q) and /3 0j i are the coefficients of 
the QCD beta function in the normalization where f3 = (11 — 2Np/3)/4. The expressions 



for r o \ and 7 9 were given in Eq. (50). This form of Uqi can be obtained from the results 
in Ref. 
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VI. COMPARISON WITH THE COLLINS-SOPER-STERMAN APPROACH 



We now compare with the standard CSS approach to transverse momentum resumma- 
tion, and demonstrate that the logarithms resumed are equivalent through next-to-leading 
logarithmic accuracy (NLL), i.e., through next-to-leading order in the resummed exponent. 
We outline what further calculations are needed to extend this result to higher orders. 

We begin by comparing the exponents that implement the resummation of large loga- 
rithms of the scales M and pr- The CSS approach writes the transverse momentum distri- 
bution as 



d 2 a 
dpxdY 



cr 



x exp 



d 2 b A 

W) 



j-ipT-b 



^ [C a ® fa/p] (%a, b /b±) [C b g> f b/ p] (x B , b /h 




a.b 
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In— -A(a s (fj,)) + B(a s (fi)) 



F 



(61) 



where we have neglected the remainder term Y which also appears for the purpose of this 
discussion. The coefficients A, B, and C a have perturbative expansions in the strong coupling 
constant: 






c„ 



E 

?1=0 



a s (jj) 



TV 



Ci n) . (62) 



The results for these coefficients are well-known in the literature [7]. 

The exponentiation of low pt logarithms in our approach is accomplished via the RG 
evolution of the hard function H from \xq ~ M to jjlt ~ Pt- As noted before in Section III A 



the matching coefficient here is the same as that for inclusive production of vector bosons 
or the Higgs. Therefore, we can take the solution for the hard-function evolution from the 
literature |31i|: 



H(Q 2 ,fi Q ,fi T ) = H(Q 2 ,fi Q ) 



""Q d[i 2 
x exp<; / — 



Q 2 

hi—r H (a 8 ([j,)) + 7i/[a s (/i)] 



(63) 



where we have set Q 2 = X\X2Q 2 - The expressions for the cusp anomalous dimension Y h and 
7 y for the cases of Higgs and electroweak gauge boson production are given in Refs. [3H [36] . 
They have perturbative expansions, as do the A, B, and C a coefficients of the CSS approach. 
A detailed study of the anomalous dimensions in SCET and in the standard QCD approach 
was previously studied in the context of threshold resummation [31]. The cusp anomalous 
dimension is equivalent to the A factor appearing in the exponent to the order we are 
working, r#[a s (//)] = A(a s (/j,)). The leading terms in the anomalous dimensions that 
controls the single logarithm are the same: ^ vl ~ l "> = B^. At the two-loop level, this is no 
longer true; the effective theory organizes terms differently than the standard approach, and 
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contributions from the matching coefficients H(Q 2 , jjlq) and Q 1 ^: j v ^ = B^+ contributions 
from H(Q 2 , /j,q), Q % K This has been observed in previous analyses comparing SCET evolution 
to the QCD literature [HI [31]. A two-loop computation of Q % i is required to further check the 
relation between the CSS and effective theory approaches; this calculation is an important 
goal for future work. By design, both approaches fully reproduce the low-pr limit of fixed- 
order result upon expansion in as- A check of the NLO pt spectrum would also require a 
two- loop computation of Q qrs . 

However, note that in the SCET approach, the low scale endpoint of the RG evolution 
of the Sudakov factor is at //y ~ px- This differs from the standard approach where the 
corresponding endpoint is at fi ~ 1/&_l where b± is the impact parameter that is integrated 
over from zero to infinity. The limit of bj_ — > oo gives rise to a Landau pole that must be 
dealt with by introducing an external prescription for any value of pt- The SCET approach 
avoids this issue as the RG evolution is done entirely in momentum space. 

A well-known aspect of the CSS approach is its treatment of the limit pt — > 0, M — > oo. 
It predicts that in this limit da jdp\ goes like a power [3] of Kqcd/Q arid is thus sensitive to 
non-perturbative input. In the effective field theory approach this corresponds to the region 
where the TMF is no longer perturbative. The leading \jp\ term coming from perturbative 
soft and collinear gluons is strongly Sudakov-suppressed by the evolution due to the cusp 
anomalous dimension in this limit. The remaining contribution then comes from the non- 
perturbative region in the effective theory whose analysis remains to be done. 

We also note that in our formalism, the factorization theorem is in terms of iBFs and 
the iSF which differ from the corresponding objects in the TMD factorization formalism 
[STH^o] . The iBFs and iSF depend on additional light-cone residual momentum components 
which regulate rapidity divergences as dictated by the physical kinematics of the process 
at finite px- This allows one to compute in perturbation theory the iBFs and iSF in pure 
dimensional regularization without the need for additional external regulators as in the TMD 
factorization approach. One can obtain a factorization formula in SCET analogous to the 
TMD factorization by expanding in these residual light cone momenta to get JT7] 

d 2 a ^ n F KL * f 2 f d 2 b ± ib ± .k ± ,\ ?2 ,,21 zrKLtijq, s 

1^ 4Q4/V2 / d k± / Jo^ 6 ± ^[^uUJt-k ± -M z ^ H z ' jq {u u ,u t ,[iQ;nT) 
x J*(u u ,0, b±,ii T )Jn( u t, 0,bx,/J 1 T)S qq (0,0,b±,n T ), (64) 



du dt 

qijKL 



where u u = Z Q U ,u t = z ~ . The SCET objects J%, J|, and S qq contain spurious rapid- 
ity divergences that require additional regulators beyond dimensional regularization. For 
perturbative values of pt, a matching calculation can be performed to write the above for- 
mula in terms of standard PDFs. A more detailed comparison of our approach to the TMD 
factorization formalism is left for future work. 
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A. Expansion of resummed formula to higher order 

To demonstrate explicitly that our formalism correctly obtains the large logarithms of the 
CSS approach at higher orders, we expand the resummed Z-boson differential cross section 



of Eq. (59) to 0(a 2 ). Our derivation closely follows the approach of Ref. [H]. We begin by 



inserting the partonic cross section of Eq. (59) into the hadronic convolution of Eq. (56) 



The resulting expression has the schematic form 
d 2 * z 



dp 2 T dY 



/ dx\d%2 F(xi, X2) 6 (x\X2-s + x\{t — M§) + X2(u — m 2 z ) + M z ) 



(65) 



where s, t, and u are the usual hadronic Mandelstam variables, defined for completeness in 
Appendix Kj The function F(x±, X2) denotes the contributions from the matrix elements and 
parton distribution functions. The delta function comes from the partonic differential cross 
section, and allows one of the integrals over partonic momentum fractions to be immediately 



performed. It is convenient to divide the integration in Eq. (65) into two regions: one where 
the pseudorapidity of the parton recoiling against the Z to give it a transverse momentum 
is greater than the Z rapidity Y, and one where it is less than Y. Doing so, and using the 
delta function to perform the X2 integration in the first region and the x\ integration in the 
second piece, leads to the expression 

d2az - r dx F{xi ^ + r dx f{x ^ x2) reel 



dp^dY J \n^e Y x i s + u — M% J /5q: e -y x 2 s + t — M 



where 
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x 2 s + t- M 2 



z 
Xl (M 2 -t)-M 2 

X\S + u — M 



2 



= h±±Ml + M (67) 

V s V s 

To proceed, we now simplify the partonic cross section by expanding around the pt — > 
limit and keeping only the l/p% terms. For simplicity we focus henceforth only on the qq 
partonic channel; the qg channel proceeds identically. In the first region of the integration, 
the partonic Mandelstam variables simplify in the px — > limit as follows: 

i->Ml(l- — ), u^O, s^ Ml—, (68) 

V x aJ x a 

where we have introduced the notation 

M Z y M Z _y 

xa = —j=z , %b = —j=z ■ (69) 

V s V s 
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The function F appearing in the integrand takes the following form in the first region after 
this simplification: 



F q q(x U X* 2 ) -> f q /p{x l )f q /p(x* 2 ) X - T 

Pt 



1 + 



Xa 
X\ 



(70) 



where for simplicity of presentation we have suppressed the overall constants which appear. 
A similar simplification and structure are obtained in the other part of the integration. 

We reduce this further by simplifying the remaining integrals over the Xi, following the 
procedure outlined in Ref. [IT]. To facilitate comparison with results in the literature we 
introduce the standard notation for the convolution of two functions, 



{f®g)(z) 



dxdyf(x)g(y)5(z-xy), 



(71) 



and remind the reader of the leading-order DGLAP kernel 

'1 + x 2 ' 



P m (x) 



C h 



1 — X 



(72) 



We also introduce the following combinations of coupling constants to match the notation 
in Ref. [IB], with which we eventually compare: 



qq 16cos 2 ^ 



[l + (l-4|e g |sin 2 ^) 2 ] 



(73) 



For simplicity we continue to focus on the qq partonic channel. After straightforward ma- 
nipulations we arrive at our result for the differential distribution: 
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, +31n^f 
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(74) 



We have explicitly denoted the scales which appear in the overall coupling constant and in 
the PDFs. We note that the solution for the evolution factor XJui can be obtained from 
Ref. [31]; to the order in as we are working, the different momentum scales which appear 
in the evolution factors in the partonic cross section do not matter, and a simple overall 
exponential factor is obtained in the differential cross section. 

To compare the structure of logarithms with those obtained in the CSS approach, we 
first use renormalization-group arguments to evolve all coupling constants which appear to 
an arbitrary renormalization scale ftp, and similarly use DGLAP to evolve all PDFs to the 



factorization scale \ip. We organize our result following the notation of Ref. 
expansion in a s and \d.(Mz/pt)'- 

d 2 a Z ,oa 4:TC 2 
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into a joint 



m 



a 



dp\dY 



3 sin 2 ^ 



</</ 



s Pt 



2tt 



Rj 



D ln-M 
Pt 



(75) 



We set /xq = Mz and /iy = pt (we comment later on the choice of an imaginary matching 
scale fiQ, as suggested recently [49]). Only terms through 0(a 2 s ) are kept. We introduce the 



explicit forms for the first few coefficients appearing in the CSS expansion of Eq. (62 ): A^ 



2Cf, B^ 1 ' = —3Cf- Introducing the nomenclature f q /p(x A , Hf) = f A , fq/p{%B, /Af) 
we find the following results for the first few coefficients: 



./' 



B, 



2 D 2 



A {1) f A f B , 

B {1) f A f B + h (P gq ® f) A + h (P qq ® f) B , 

[AU] 2 f A f B , 
\AV> [f B (P m ® f) A + U (P„ ® f) B ] - 



1 

2 

3 

2' 



{-AW/b (P gg <g> /)a^ " 25(1) /* tf« ® /)* " ^ [ 5(1) ] 2 /a/b 



A«B« - f3 A^ 



/a/, 



Bj 



A 



A> 



+ |^) /A/B lng| + ^)/^-(P f . 
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/)A (^ ® /) 



B 



-f B (P qq <g> P OT ® /) A + /3 /b (P m ®f) A } + [A^B}. 



(76) 



The coefficients i-Di, iD , 2 B>3, and 2-D2 agree 2 with the analogous n C m coefficients of 
Ref. jH] that appear in both the fixed-order expansion and the CSS formalism. Differ- 
ences occur in 2 -Di; the 2 C\ formalism of the usual approach contains two additional terms 
depending on the quantities A^ and C^ 1 '. This is not surprising, as our result has been 
computed only to next-to-leading logarithmic accuracy. These terms in the expansion are 
of next-to-next-to-leading logarithmic order. Denoting L = In M\jp\^ we remind the reader 
that resummation to a given order gives the following towers of logarithms [10]: 



leading logarithmic 

next-to-leading logarithmic 

next-to- next-to-leading logarithmic 



^ n T 



n r 2n— 1 

a?L 2 "- 2 , 



n j 2n— 3 



^ n t 



(77) 



The full result at next-to-next-leading logarithmic accuracy along with the complete result 
for 2-D1 requires the next higher order calculation of the TMF. However, some of the next-to- 
next-to-leading order logarithmic terms can be already seen to appear in the partial result 



2 We disagree with the statement made in Ref. |50j that our formalism does not correctly resum logarithms 
at the next-to-leading- logarithmic order; our explicit check makes it clear that this claim is incorrect. 
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for 2-D1 above. In our factorization formula, we obtain two distinct large logarithms: explicit 
logarithms ln/ig//x^ coming from the resummed exponent, and a kinematic one of the form 
lnM|/py coming from the hadronic convolution. This organization is different than in the 
CSS approach, but all required terms are correctly obtained. 

It was recently suggested in the literature to utilize an imaginary matching scale for fig, 
which has the effect of resumming factors of n 2 which arise from the time-like momentum 
transfer [36, 49J. This was shown to improve the convergence of the perturbative expansion 
for inclusive Higgs production [361 EH] , and has also been utilized in the literature to study 
Drell-Yan [35]. Doing so here adds the following additional term at the single-logarithmic 
order: 

2 D 1 ^ 2 D 1 + C F TT 2 f A f B . (78) 

This factor is part of the contribution to the coefficient C^- 1 ' in the CSS approach. 

VII. NUMERICAL RESULTS 

We present here numerical predictions utilizing the factorization and resummation for- 
mulae we have derived. We show results for Higgs production at a 7 TeV LHC, and for Z 
production at the Tevatron. Our numerical results are based on the resummed partonic cross 



sections presented in Eq. (59). For the Z we compare with the Run 1 data from both CDF 
and DO to demonstrate the consistency of our calculation with experimental results. Our 
results are model independent and free of Landau pole prescriptions required in the standard 
approach. For perturbative values of pt, the transverse momentum distribution is given en- 
tirely in terms of field-theoretically derived perturbative functions and the standard initial 
state PDFs. For non-perturbative values of p?, the TMF function Q is non-perturbative, 
but field-theoretically well-defined, and can be modeled and extracted from data. In this 
section, we restrict our results only to perturbative values of pt leaving the non-perturbative 
region for future work. 

Before describing the parameter values assumed in our study, we comment on the order 
of our resummation. The hard matching coefficient H, the cusp and non-cusp anomalous 
dimensions Th and 7# respectively, and the TMF Q all have perturbative expansions that 
must be calculated to sufficiently high order to achieve resummation of certain classes of 
logarithms. We note that since generating a finite pt requires the emission of at least one 
parton, the contribution of the TMF to the transverse momentum spectrum begins only 
at 1-loop. All quantities required to achieve next-to-leading logarithmic accuracy (NLL) 
are known, and have been detailed in previous sections of this paper. To achieve NNLL 
precision, the 2-loop result for Q is needed. 

We start by discussing some general features of our numerical results. We show fixed-order 
results at leading order in perturbation theory, and results at both LL and NLL matched 
to the fixed-order results at (D(a s ), as shown in Section [Vl In the standard nomenclature 
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VS = 7 TeV 
p T 2 /2 ^ fx T z < 2p T 2 
-m h 2 /2 < Mq 2 ^ -2m h 2 
= 165 GeV 




50 75 100 

Pt (GeV) 



150 



FIG. 4: Numerical predictions for the transverse momentum spectrum for Higgs boson production 
at the LHC for central rapidity. Shown are the fixed-order result and those obtained after imple- 
menting the resummation formula of Eq. (J™ through LL and NLL. The bands arise from the scale 
variation shown in the text. 



these would be called LL+LO and NLL+LO. We use MSTW 2008 parton distribution 
functions [5T] . For LL and LO predictions we use leading order PDFs with 1-loop running 
of the strong coupling constant, while for our NLL results we use NLO PDFs with 2-loop 
running for a s . Our results depend on the two matching scales [A? and /iq. The dependence 
of the cross section on these arbitrary scales occurs at one order beyond the order in a s to 
which we have calculated; it would vanish completely if we could compute the cross section 
to all orders in perturbation theory. The variation of these scales therefore provides some 
indication of missing higher-order effects, and is conventionally used as an estimate of the 
theoretical uncertainty. We must choose both a central value for these scales and a range 
of variation to obtain an uncertainty estimate. As our central scale choices we set /i^ = p\ 
and Hq = —M 2 . These are chosen to minimize logarithms that appear in the perturbative 
expansions of the hard function and the TMF, as discussed in Section pi] We vary /z^, 
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FIG. 5: Numerical predictions for the transverse momentum spectrum for Z boson production 
at Tevatron Run 1, compared with data form both CDF and DO. Shown is the resummation 
prediction of Eq. (pi) at NLL. The bands arise from the scale variation shown in the text, while 
the result for the central scale choice is shown by the solid line. The lower limit of the plot is pt= 
1.75 GeV. 



/Xg independently around these choices by a factor of 2. Two unconventional aspects of 



these choices require comment. Following Ref. [3S], we utilize an imaginary matching scale 
for fiQ, which has the effect of resumming factors of rr 2 which arise from the time-like 
momentum transfer appearing in H . This was shown to improve the convergence of the 
perturbative expansion for inclusive Higgs production (36J, 09], and has also been utilized 
in the literature to study Drell-Yan [35J. We also find better agreement with data (see 
Fig. [5]) for an imaginary fig compared to a real jjlq which can be attributed to the effect of 
resumming factors of tt 2 with the former choice. We also choose to vary our scales around 
a reduced range to avoid evaluating a s (nT) at a non-perturbative scale when the transverse 
momentum becomes small. An framework for incorporating the non-perturbative region of 
transverse momentum into the SCET formalism was given in Ref. [52] • In this approach the 
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scale (jlt freezes at a value near Aqcd as the px approaches zero. 

In Fig. [4] we show the predictions for the Higgs px spectrum at the LHC, using both 
the fixed-order expression and the resummed results at LL and NLL accuracies. The 
general features of this plot are clear: large logarithms of the form In (m\/p T ) spoil the 
fixed-order perturbative expansion at low px- The Sudakov suppression coming from the 
renormalization-group evolution of the hard function H tames this behavior. The central 
value of the prediction is absolutely stable upon proceeding from LL to NLL; only a reduc- 
tion of the scale variation is observed. At intermediate and high momenta, the matching 
onto the fixed-order expression is smooth. The sensitivity to scale choices that can lead to 
negative results [15] in the standard approach, does not occur in this effective field theory 
approach. This allows the matching scale jiq to be varied throughout a range sufficient to 
use it as an estimator of the theoretical uncertainty. An additional error also arises from 
imprecise knowledge of parton distribution functions. We postpone a numerical analysis of 
this issue until a detailed study of boson p T distributions at both the Tevatron and LHC 
incorporating the non-perturbative region is performed. 

One aspect of transverse resummation in SCET that requires further study is the treat- 
ment of the non-perturbative region px ~ Aqcd- in our analysis, the transverse momentum 
function Q becomes non-perturbative, and must be modeled. The onset of this region can 
be seen in the plot by the large scale variation at low px, which is caused by evaluating 
&s{ht) ~ o: s (Aqcd)- Since this object has a non-pertubative definition in terms of operator 
matrix elements and a well-defined running, it is reasonable to extract this function using 
available data. In our plot for the Higgs px distribution we simply stop our plot at a lower 
value of px = 3 GeV. The study of the non-perturbative region of px was recently begun in 
Ref. [52]. 

In Fig. [5] we plot our prediction for the Z-boson px distribution at the Tevatron Run 1, 
and compare to data from CDF [21] and DO [23] • We study the spectrum down to px = 1-75 
GeV. The agreement with the data is excellent over the entire range. The low px version 
of this data can eventually be used to constrain the non-perturbative TMF that appears in 
SCET, as is done in the CSS approach [53] . 

VIII. CONCLUSIONS 

In this manuscript we have extended our analysis of transverse momentum distributions 
using the Soft-Collinear Effective Theory(SCET) to account for both electroweak and Higgs 
boson production at low px in hadronic collisions. We have derived a factorization theorem 
for the transverse momentum distribution for the production of electroweak gauge boson 
production, and have provided all necessary analytic expressions to perform resummation 
of \ow-px distributions for any color-neutral particles to next-to-leading-logarithmic accu- 
racy. Our effective field theory approach is free of the Landau pole that appears in the 
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standard approach even for perturbative values of px- We thus have a numerically stable 
matching to the fixed-order QCD result, leading to a smooth transition from the low-pr 
resummation region to the intermediate and high p? region without the need for a matching 
prescription. For perturbative values oipT, our approach predicts the transverse momentum 
distribution entirely in terms of field-theoretically derived perturbative functions and stan- 
dard initial state PDFs. For non-perturbative values of px, an additional non-perturbative 
Transverse Momentum Function (TMF) appears with a rigorous field-theoretic definition 
and computable anomalous dimension. 

We have presented the first numerical predictions for pt spectra arising from SCET for 
Higgs and Z-boson production, and for Z boson production have shown an initial numerical 
comparison with Tevatron data. The agreement with data is excellent over the kinematic 
range currently covered by our factorization formula, indicating that SCET will provide a 
useful framework for the analysis and interpretation of hadron collider distributions. 

Our analysis reveals several directions for future work. Precision predictions at next-to- 
next-to-leading logarithmic accuracy require two-loop computations of the iBFs and iSFs 
that appear in our factorization theorems. This computation is within current technical 
capabilities. The region of non-perturbative pt requires further study through a modeling 
of the non-perturbative TMF followed by its extraction from data. 

In summary, the SCET approach to transverse momentum resummation offers a com- 
pelling alternative to the usual method. Several explicit checks have been performed: (1) 
a comparison with the leading fixed order cross-section, (2) the cancellation of the scale 
dependence between the various components in the factorization theorem as determined by 
their RG evolution structure, and (3) an explicit check of the logarithms at next-to-leading 
logarithmic order. Furthermore, we find good agreement with data. We look forward to the 
further development of our results. 
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Appendix A: Factorization of electroweak gauge boson differential distributions 

In this appendix we describe the steps in the derivation of the factorization formula for 
the transverse momentum and rapidity distributions for electroweak gauge boson produc- 
tion. These steps closely follow the derivation of the analogous factorization formula for 
Higgs production derived in |19| . For simplicity in notation we focus here on the case of 
single-boson production. With straightforward modifications, one can obtain analogous fac- 
torization formulae for neutral-current 7* + Z production as well as the case where the final 
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state leptonic decay products of the vector boson are treated differentially. It is convenient 
to work with the hadronic Mandelstam variables u, t which are related to pt and Y as 



u=(p2- q) =M 2 Z - QJp 2 T + M 2 e 



t = (pi- qf = Mf - Q^Jpl + M 2 e~ Y , 

(Al) 

where q^ and M z denote the vector-boson momentum and mass respectively and 3 

dudt = Q 2 dp 2 T dY. (A2) 

After matching the vector-boson production current onto SCET PT current as explained 



in section III A , the differential cross-section in the hadronic Mandelstam variables takes the 
form 



d 2 a I rl-i f d 4 q 



i27i)6(q°) S(q 2 - M 2 z )L» v {q) I du x I du 2 f dcu[ I dJ 2 



dudt 2Q 2 L4J J (2vr) 4 

x E E (^) 4 ^ ) (Pi+P^-q-Px n -Px fi -Px s )S[u-(p 2 -q) 2 ]S[t-(p 1 - q ) 2 ] 

qq'ijKLX„,Xn,X s 

x C'^( Wl , W2 )C'* L ^\^,^)(p 1 p 2 |O^'tK,4)(0)|X n ,X n ,X s ) 

(A3) 

where the indices run over 

K,L = {V,A}, q, q ',i,j = {u,d,s,---}, (A4) 

where V and A label the vector and axial-vector Dirac structure. The indices q,q',i,j run 
over the massless quarks that appear in the initial protons. The contribution from a pure 
gluon SCET PT operator that would produce gg — )■ V with V = 7*, Z, W vanishes for Drell- 
Yan processes, so that the sum over i,j does not include the gluon. The overall factor of 
I comes from averaging over the initial hadron spins, the final state \X) has been broken 
up into the n-collinear, n-collinear, and soft states so that \X) = \X n X n X s ), and the SCET 
operators have the form 

(A5) 
where the Dirac structures T^ q are given by 

r v , q = ghi, r? = 9hiis- (A6) 

3 Note that Q = n ■ p\ — n ■ P2 denotes the hadronic center of mass energy, and is not related to g M which 
is the vector-boson momentum and satisfies q 2 = M%. 
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g v and g\ denote the vector and axial- vector couplings of the q-th quark to the vector boson 
respectively. The tensor L^ v denotes the product of the leptonic currents arising from the 
vector-boson decay. For simplicity of notation, we will present out formulae integrated over 
the leptonic phase space, so that we can use effective polarization vectors and set 



V»(q) = J2 e/ W(?) = -gT + q -X- (A7) 



M: 
pols. 



By the equations of motion for massless quarks, the contribution of the q^q v term in the 
above polarization sum vanishes when contracted with the quark bilinear currents. This 
allows us to effectively set L v u (q) — y —g^ u which is used in the rest of this analysis. Next we 
use the soft-collinear decoupling |2D1 EI] property of the leading order SCET PT Lagrangian 
to decouple the matrix elements into n-collinear, n-collinear, and soft objects so that the 
differential cross-section becomes 



(Pa -g^ f d A q 



%°) 5(q 2 - M\)b [u - (p 2 - qf] 5[t - ( Pl - qf] J du 2 J dw[ f dJ 2 



dudt 8Q 2 J (2tt) 3 ' 

x E E ^Y^\pi+P2-q-Px n -Px,-Px s ) 

qq'ijKLX n ,Xn,X a 

x C**(c*,o*)C7^ViV2)(I?*W 

(A8) 

We perform a series of steps that allow us perform the sum over the states X n , X n , X s while 
consistently maintaining the final state restriction on the gauge-boson momentum. We begin 
by inserting the identity operator 

1 = J 'd* Pn j ' d'p n j ' d%5^\ Pn - P Xn )8^(p n - p Xn )5^{ Ps - P Xa ), (A9) 

and decompose the momenta into label and residual parts so that 

X n X n X n ' X n Xfi Xn ' V ' 

where P x , P x _ ~ Mz and K x , K x _ <C Mz- We write the remaining momentum compo- 
nents as 

Pxn=^~xt^ P xt = K 'x^ Px s =K Xs (All) 

and similarly write 



Pn — Pn + Si ' Pn ~ Pn + % ' 

Pn = k n ' , Pn' = k n ' p s = k s , (A12) 
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where again p~,p^ ~ Mz and k~,k£ <C Mz- The delta functions in Eq. (A9) can be broken 
up into Kronecker deltas over label momenta and residual delta functions which we can 
write using the integral representation as 

1 = J2 6 p-.p- ¥ P + I d A k n d A k n d A k s 6 (4 \k n - K Xn )S^(kn - K Xn )5^\k s - K Xa ) 

Pn,Pn 

/f d A x f d 4 v f d 4 z 
d A k n d A k n d A k s / — —- / — — - / — — - 
J (2vr) 4 7 (27r) 4 i (2tt) 4 



i(k n -K Xn )-x i(k n -K Xn )-y i{k a -K Xa )-z 



(A13) 



Similarly, the momentum of the vector boson can be divided into label and residual compo- 
nents so that 



n-q = n-q + n-k, n-q = n-q + n-k, q± = kj 



(A14) 



where the n ■ q,n ■ q ~ Mz and n ■ k,n ■ k,k± <C Mz- The phase space integral over g M can 
now we written as 

d A q5(q 2 -M 2 z ) = J2 f d 2 k ± f d ^^5{q + q- +q + k- + q-k + + k + k- -kl-M 2 z ). 



<? T ,9" 



(A15) 



The four-momentum conserving delta function appearing in the differential cross-section can 
be written as a product of Kronecker delta functions over label momenta and delta functions 
over residual momenta as 

5 {4 \pi +P2-q-Px n - Px* - Px s ) = 8 Wl4 -5^ A+ 8^\K Xa + K Xn + K Xn + k ± ) 

x S(K Xn + K Xn + K Xa + k + )5(K Xn + K Xn + K Xs + k~) 

= WW^W + k^ + ki + k ± ) 

x 8(k+ + ki + kf + k + )5(k~ + k^ + k; + k~), 

(A16) 



where we used the residual delta functions in the first line of Eq.(A13) to obtain the second 



equality above. Using Eqs. (A13), (A15), and (A16) the cross-section can be brought into 
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the form 

£a - =^l E E / <^ / ^%-v + *+*- + r *+ + *+*- - Mi) 

qq'ijKL q+ ,q~- 

x e i ^ x e^e*'-** WM --<W ^(^ + ^ + ^ + fc ± ) 
x 5(A;+ + A;+ + fc+ + fc + )<J(A;- + A£ + fc; + h~) 

x VWrf% ff (r^'t)^< Pl |(^^ 

x Tr(0\T[(Siy%S n y\z)}T[(Styf(S n y d }\0}, (A17) 

where we have simplified the color structure using the identities 

Xqq'fiad 

{pM q w)T{x)\x n ){x n \{wH q y n iM\p^) = - 1 ^(pimwyJ(x)\x n )(x n \(w^ l y n i i (o)\pi) 

Xqq' Xbe 

(P2\(W^ q ,t\y)\X n )(X n \(Z q W)- 2 (0)\P2) = ^{P2\(W*Q? f (y)\X n )(XMW)tM\P2)- 

(A18) 

Next we apply a spin Fierz identity which allows us to bring the cross-section into the form 
d 2 a v^ nF KL ^ v-^ f , 2j [ dh + dk~ ,,„,„ „+, „ I+ 7+7 ,^ 

dwi / du 2 H z L ' m ( y u 1 ,u 2 , fig] fi T ) / d^knd^knd^ks j -— j ——^ j ——^ 

x e i ^V^V fc '-*<L 1)9 --<W< j(2 W + ^n + fcfl + ^)5(A;+ + A;+ + fc+ + fc+) 

x «J(fc- + A;J + fcj + k-)5[u - M z + Q n- q]S[t - M z + Q n- q] 

x J q (u u x, /j, T )J?(u 2 ,y, /j, T )S qq (z, /j, t ), (A19) 

where we have defined the hard function 

HZ L * jq (u u u 2 ,ii) = C k ™(u u u 2 ,ii)C* l ™(u u u 2 ,ii), (A20) 

and H z ' ljq (oji,oj 2 , \Xq\ Ht) denotes the RG-evolved hard function from the scale \xq ~ Mz 
to [It ~ Pt- The quantity F KL ' ,t comes form the contraction of the leptonic tensor with the 
Dirac structure of the hadronic tensor, and is given by 

pKL;q _ au il__fl_l__v__i (A21) 



33 



The jet and soft functions are given by 

J>,X,// T ) = J2 (p^ W U^p{Pn-0o)(W^ q ) n {0)\ Pl }, 

initial pols. 

4(u 2 ,y,fi T ) = J2 (P^s P in[^(W^Q n (y)5(-co-P n m g WU0)]\p 2 ), 

initial pols. 

S qq (z,pL T ) = Tr(0\f[StSn}(z)T[SlS n }(0)\0). 



(A22) 



Next we perform the integrals over x + and y components in Eq. (A19) by defining the 
Fourier transformed jet functions as 

dx ife— x + 



dy~ 

Air 



ea nX J%(u 1} x + ,x ,x ±)/ u) = J^(ui,k n ,x ,x±,n), 



-ea 



Kv Jn( UJ 2,y + ,y ,y±,ti) = Jn( UJ 2,y + ,kf,y ±J fi) J 



(A23) 



combining the label and residual momenta as 



doj\dk n — ¥ / dui, 



dtu 2 dk£ — > / du 2 



(A24) 



and absorbing the residual momenta k n , k£ into u\,u 2 respectively to get 



(Pa 



E 



7T 



F KL ' C1 f dq+dq- 



dudt ^ 4Q 2 A" r 2 / 2 

qijKL ^ c 



d 2 k ± 5(q + q -k\-M\) \ du>i / du 2 H z '' l39 (u 1 , co 2 , I^q] Ht) 



dx d 2 x± f dy + d 2 y± 1 f dz + dz d 2 z± f db + db f d 2 b± 1 



(2tt) 3 J (2vr) 3 2 J (2vr) 



2(2tt)2 J (2tt) 2 2 ' dft ™ dft ™ 



(A25) 



x \ I dk^d 2 k± l - f dkjdk;d 2 kje- i ^<^-^ ) e- % ^-^-^ ) e- t ^<^-^ ) 

x e ifc+(x--b-) e |fcr(y+-6+) e ifc7(^+-b+) e ffc+( z --fe-) e i( Wl - g -)6+ e f( W 2-g + )b- e ^±-fc± 

x 5 [m - Mf + Q g _ ] <J[t - M 2 Z + Q g + ] J*(wi, x", x ± , /i T ) J|(w 2 , y + , y±,n T )S qq (z } /i T ). 

Performing the integrals over the momenta k^,k^,k^ n and fc^ and the x,y,z coordinates, 
we get 

S5 = £ dW ^1^ %V - £ - «S)*[« - Mi + Q r]t[t - Ml + ,+] 



2(2tt) 2 J (27^ 
x J^(wi, b~,b ± ,n T )J?(uj 2 , b + , b ± ,fx T )S qq (b + , b~,b ± ,fx T 



(A26) 
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where from brevity we have defined 



flfKwi.AkjJfir) = J2 FKL ' qH z L]lJq ^i^2^ Q ^T). (A27) 

KLij 



This expression can be brought into the form 
cPa *-^ vr f , , , /" , Q , /" cPb 



du dt 



^ Wm I dq+dq ~ I ^ I J^T^ ±rk±6[u -M 2 z + Q q -]S[t-M*+Q g+] 
x £(<? +, T - fc± - M|) / dux / dw 2 (4wiu; 2 ) H q z (ui,u 2 ,iJ,Q\ Hr) 
x / dk*dk^Bl{u Xl fc+, 6 ± , fi T )Bl(u 2 , k^, b ± , ti T )S qq {uJi - q~ - k^,u 2 - q + - fc+, 6±, /x T ) 



(A28) 



where we have made use of the Fourier transformed functions defined as 
ASK #,&.!.,/*) = ^j-J^e^-J^b-^ri, 
Bl(u 2 ,k^,b ± ,n) = — — e5^ fe+ j|( W2 ,6 + ,6 ± ,/i), 
5 w (wi,w 2> 6 J . )A i) = f d ^ r e^ b+ e^ h -S qq {b + X^^). (A29) 



The -B,^'! functions are referred to as the purely collinear impact-parameter Beam Functions 
(iBFs) and are defined with a zero-bin subtraction to avoid double counting soft emissions 
already encoded in the soft function S qq . It was shown in [TJI] that convolution over the 
purely collinear iBFs and the soft function can be written as 

x <$(<? + <r ~k\- M z ) j dux j du 2 (4wiu; 2 )#| (wi,u; 2 , /jLq; /Mr) 

x / dk^dk^B^(uj 1 ,k^,b ± )Bl(uj 2 ,k^,b ± )S qq l (uj 1 - q~ - k n ,u 2 - q + - k+,b ± ), 

(A30) 

where B% q n are the 'naive' iBFs or simply the iBFs defined without a soft zero-bin subtrac- 
tion, and S^ 1 is the inverse Soft Function (iSF). Next we rewrite the cross section in terms 
the variables xi,x 2 ,t+,t^ defined as 

xi = -^, x 2 = -^, i+ = xiQfc+, tr = x 2 Qk^, (A31) 
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to get 
d 2 a 



71 



dq dq Id k 



d 2 b 



1 Jbj -k 



e*^' fc± 5[u - M 2 Z + Q q-]5[t - M 2 Z + Q q + ] 



dudt N 2 J * * J "J (2tt) 2 

x 5{q + q~ - k\ - M 2 Z ) / dx\ / dx 2 V^ H q z (x 1 x 2 Q 2 , Hq] /xt) 

Jo Jo 

x / dt+dt^B%(x u t+, b±)Bl(x 2 , t^, b±)S~ {xiQ - q~ ^,x 2 Q - q + ^, b 



x 2 Q' 



xiQ' 



(A32) 



where we used the fact that H z {uji,uj 2 , Hq', />r) — H z (cuiu>2, Hq', />r)- in the next step, the 
iBFs are matched onto the PDFs as 



Bn,n{x,t,b ± ,fi) 



■J-n,n;qr\ i^il^) Jr\Zj f*)i 



(A33) 



so that the differential cross-section becomes 



d 2 a 



TX 



du dt Q i\ c J j j Xl x\ j X2 x' 2 



dxi I dx 2 



(jjJU 1 / (Jitlj Q 



X 



^2 H z (xiX 2 Q 2 , fi) G qrS (x!, X 2 , x[, X 2 , U, t, ^ T )fr{.x' 1 ^ T )fs{.x' 2 , /i T ), 



(A34) 



with a sum over repeated indices understood. The function Q qrs is given by 

d 2 b 



Q qrs (x 1 ,x 2 ,x 1 ,x 2 ,u,t^ T ) 



(M 



-Jo 



hKi { Ml-u)iMl-t) _ Ml 



X / dt^dt n l n - qr ( — ,t^,b_ L ,^ T )In;qs{ — ,tn,b±,HT 

J x 1 x 2 

(M§ - u) t, 



x Sj{x\Q 



Q 



x 2 Q 



,x 2 Q - 



M 2 -t t+ 



Q xiQ 



(A35) 



Using Eqs.(Al) and (A2) we can obtain the differential cross-section in terms of the px and 
Y variables 



d 2 a 



,2 /•! 







;r i " r r dx[ f dx 2 

ax i / dx 2 






dp 2 T dY N 2 J{) 

x 5Z H z( x i x zQ 2 i A*Q5 A*t) G grs (x u x 2 , x[, x' 2 ,p T , Y, n T )f r (x' v /j, T )f s (x 2 , /i T ), 



(A36) 
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where 

Q qrs (x 1 , x 2 ,x' 1 ,x 2 ,p T , Y, fpp) 



J (2tt) 2 J x\ x' 2 



x S~ 1 (x 1 Q - e Y Jp 2 T + Mf - ^,x 2 Q- e- Y J V \ + Mf - -^-, b ± , fi T 



x 2 Q v xiQ 
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